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Abstract 

Despite its century-old use, the interpretation of local field potentials (LFPs), the low- frequency part 
of electrical signals recorded in the brain, is still debated. In cortex the LFP appears to mainly stem 
from transmembrane neuronal currents following synaptic activation, and obvious questions regarding 
the 'locality' of the LFP are: What is the size of the signal-generating region, i.e., the spatial reach, 
around a recording contact? How far does the LFP signal extend outside a synaptically activated neu- 
ronal population? And how do the answers depend on the temporal frequency of the LFP signal? 
Experimental inquiries have given conflicting results, and we here pursue a modeling approach based 
on a well-established biophysical forward-modeling scheme incorporating detailed reconstructed neuronal 
morphologies in precise calculations of population LFPs including thousands of neurons. 

The two key factors determining frequency dependence of LFP are (1) the spatial decay of the single- 
neuron LFP contribution and (2) the translation of synaptic input correlations into correlations between 
single-neuron LFP contributions. Both factors are seen to give low-pass filtering of the LFP signal power. 
For uncorrelated input only the first factor is relevant, and here a modest reduction in the spatial reach 
is observed for higher frequencies compared to the near-DC value Hz) of about 200 /im. Much 
larger frequency-dependent effects are seen when populations of pyramidal neurons receive correlated 
and spatially asymmetric inputs (basally or apically): the low-frequency (~ Hz) LFP power can here 
be an order of magnitude or more larger than the LFP power at, say, 60 Hz. Moreover, the low-frequency 
LFP components are found to have larger spatial reach and extend further outside the active population 
than high-frequency components. Our numerical findings are backed up by an intuitive simplified model 
for the generation of population LFP. 

Author Summary 

The first recording of electrical potential from brain activity was reported already in 1875, but still the 
interpretation of the signal is debated. To take full advantage of the new generation of microelectrodes 
with hundreds or even thousands of electrode contacts, an accurate quantitative link between what is 
measured and the underlying neural circuit activity is needed. Here we address the question of how 
the observed frequency dependence of recorded local field potentials (LFPs) should be interpreted. By 
use of a well-established biophysical modeling scheme, combined with detailed reconstructed neuronal 
morphologies, we find that correlations in the synaptic inputs onto a population of pyramidal cells may 
significantly boost the low-frequency components of the generated LFP. We further find that these low- 
frequency components may be less 'local' than the high-frequency LFP components in the sense that 
(1) the size of signal-generation region of the LFP recorded at an electrode is larger and (2) that the 
LFP generated by a synaptically activated population spreads further outside the population edge due 
to volume conduction. 

1 Introduction 

The measurement of electrical potentials in the brain has a more than hundred year old history fl]. 
While the high-frequency part has been successfully used as a measure of spiking activity in a handful 
of surrounding neurons, the interpretation of the low-frequency part, the local field potential (LFP), has 
proved more difficult. Current-source density (CSD) analysis of multisite LFP recordings across well- 
organized layered neural structures such as cortex and hippocampus, was introduced in the 1950's 
However, even if the CSD is a more local measure of neural activity than the LFP |3-8 , the interpretation 
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in terms of underlying activity in neural populations is inherently ambiguous [9|[l0]. Thus in many in 
vivo applications, for the example when investigating receptive fields in sensory systems, the LFP signal 
was discarded altogether. The LFP signal has seen a revival in the last decade, however. This is due to 
the rapid development of new silicon-based microelectrodes now allowing for simultaneous recordings of 



LFP at tens or hundreds of contacts 11 14 (and availability of affordable high-capacity hard discs to 



store the data), the realization among neuroscientists that the LFP offers a unique window into neural 
activity at the population level [9 15-23 , and the possibility of using the LFP signal in brain-machine 



interfaces 24 - 27 



To take full advantage of the opportunities offered by this new recording technology, a precise un- 
derstanding of the link between the recorded LFP and the underlying neural activity is required. For 
example, two obvious questions regarding the 'locality' of the LFP that need quantitative answers are: (1) 
What is the size of the signal-generating region, i.e., spatial reach, around a recording contact? (2) How 
far does the LFP signal extend outside an active population due to volume conduction? The first question 
has been addressed in several experimental studies, with resulting estimates for the spatial reach in cor- 
tex varying from a few hundred micrometers to several millimeters ^28j-|33j . This large range in reported 
experimental estimates presumably reflects that the spatial reach depends strongly on the spatiotemporal 
properties of the underlying spiking network activity, in particular the level of correlations [34'. These 
critical network features will not only vary between the different brain regions and species studied, but 
also depend on the brain state. 

In cortex, thousands of neurons contribute to the LFP, making the signal inherently difficult to 
interpret. Fortunately, the "measurement physics", i.e., the biophysical link between neural activity and 



what is measured, is well understood: According to well-established volume-conductor theory (10 35 



the recorded LFPs stem from appropriately weighted contributions from transmembrane currents in the 



vicinity of the electrode contact. Building on pioneering work by Rail in the 1960's 35,36 , a forward- 



modeling scheme incorporating detailed reconstructed neuronal morphologies in precise calculations of 
extracellular potentials, has been established fST' and used to explore both spikes [37-41 and LFPs f9','34| 
[41|-,43j generated by single neurons 1 37 - 40 , 42^ and neural populations [9j|34j 41 . Unlike in experiments, 
this modeling scheme allows for a clear separation between volume conduction effects and effects of 
spatiotemporal variations in spiking network activity in determining population LFPs. In f34' it was used 
in a thorough investigation of the locality of LFP. It was found that the size of the LFP-generating region 
depends on the neuron morphology, the synapse distribution and correlations in synaptic activity. For 
uncorrelated activity, the LFP represents neurons in a small region (that is, a few hundred micrometers 
around the electrode contact), while in the case of correlated input the size of the generating region is 
determined by the spatial range of correlated synaptic activity and could thus be much larger. Specifically, 
it was found that correlated synaptic inputs onto either the apical or basal dendrites of a population of 
pyramidal neurons could give orders of magnitude larger LFPs, and a much larger spatial reach, compared 
to the situations with (1) the same correlated input spread homogeneously over the neuronal dendrite or 
(2) similar uncorrelated synaptic inputs placed evenly or unevenly over the neurons. 

As shown in ^34^ , the relative contributions to the population LFP from neurons at different distances 
from the electrode will depend on three factors: First, a single neuron close to the electrode will contribute 
more to the LFP than if it was placed further away. Second, for a disc-like population, characteristic 
for a laminar population in a cortical column, it follows that with constant neuron density, the number 
of neurons located on a ring at a particular radial distance r from the electrode will increase linearly 
with r. Third, with correlated synaptic inputs onto a neural population, the LFP contributions from 
different cells will also become correlated, or synchronized, and effectively boost the contributions to the 
LFP. The contributions from different rings of neurons will thus be determined by the interplay of these 
three factors. In [34* a simplified model for LFP generation based on these elements, (1) the decay of the 
single-neuron contribution with the distance from the electrode, (2) the population geometry, and (3) the 
correlation of LFP contributions from individual neural sources, was constructed. We found this simple 
model to not only give qualitative insight into the generation of population LFPs, but also quantitatively 
accurate predictions. 

Strong frequency dependencies have been observed both in the tuning properties (28 29 and infor- 
mation content fT8','22' of cortical LFPs. For example, the low-frequency LFP (less than 12 Hz) has 
been shown to carry complementary information to the gamma-range LFP (50-100 Hz) in VI of macaque 
monkeys during naturalistic visual stimulation (22]. To properly interpret such experiments, it is thus 
important to know how spatial reach of the LFP varies across frequencies and whether the biophysics of 
LFP signal generation boost some frequencies compared to others. The high-frequency LFP components 
are, for example, expected to be more local than the low-frequency components due to intrinsic dendritic 
filtering (42], i.e. due to the reduction of the (effective) current-dipoles with increasing frequency resulting 
from the capacitive properties of the dendritic membrane |10 . 

In (34] we used the biophysical forward-modeling scheme to investigate the total population LFP, 
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i.e., the total signal generated across all frequencies. Here we use the same scheme to investigate both 
the distribution of the power of synaptically generated LFP between different frequency bands and the 
frequency dependence of the locality of the LFP signal. In terms of the latter, we study the size of the 
signal-generating region (spatial reach) as well as the spatial extension of the LFP signal outside an active 
population — for each frequency component separately. 




Figure 1. Spatial reach of different frequency components of LFP for different levels of 
synaptic input correlations Ci„. Color lines denote parts of the whole population (gray, radius = 
1 mm) which contribute 95% of LFP amplitude at given frequency in the middle of the population, at 
the soma level. Results for layer-5 pyramidal cell with basal input. 



We also use a frequency-resolved version of the simplified model developed in 34 to guide our in- 
vestigation of this frequency dependence. The population geometry (factor 2) does obviously not change 
with frequency. In contrast, the single-neuron LFP contribution (factor 1) decays faster for higher LFP 



frequencies 40 42 , but an equally important factor turns out to be the frequency dependence of the 
'correlation transfer', i.e., how correlations in the synaptic input are transferred to correlations between 
the single-neuron LFP contributions (factor 3). As an example. Figure 1 illustrates how the frequency- 
resolved spatial reach varies with the input correlation for a pyramidal population receiving basal synaptic 
inputs. We show that when the frequency dependencies of factors 1 and 3 are incorporated, the simplified 
model can still account well for the results obtained by comprehensive numerical investigations. To allow 
for direct use of the simplified model in future applications, we here thus present and tabulate numerical 
results for the frequency dependence of these key factors for a variety of situations 

The paper is organized as follows: first we describe our simulation setup, present the simplified model 
of the population LFP, and review its ingredients. Then we present detailed results of the simulations: 
we analyze the frequency content of the population LFP, the reach of different frequency components, 
and the decay of the signal outside of the population. Next we discuss the implications of our results for 
interpretation of electrophysiological data in terms of the underlying understanding the neural activity. 
Finally, in Methods we give details of the simulation setup and the mathematical model. 



2 Results 



2.1 Simulations 



The fundamental formula relating neural activity to the generation of the extracellular electrical potential, 
including the LFP signal, is given by (l0l|37| 



<^(r,i) = ^^E^, (1) 

47rcrcond ^ k - I"j 

where denotes the transmembrane current in a neural compartment i positioned at rj, and the ex- 
tracellular conductivity is denoted by CTcond- The transmembrane currents are calculated by means of 
standard multicompartmental modeling techniques with the simulation tool NEURON [U . 

An essential part of the present work is the numerical simulation of the LFP in the center of a disc-like 
population of cortical cells. The simulation setup is illustrated in Figure [2j We consider a population 
of A^=10000 cells distributed homogeneously on a planar disc with a radius of 1000 /im. Figure [2|\. The 
somas of the cells are positioned at the same depth, and the LFP is calculated at the soma level. In this 
setup we investigate how the LFP signal increases as contributions from more and more distant neurons 
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Figure 2. Simulation setup. A. Model cells are placed with constant area density on a disc of radius 
R, with the recording electrode at the population center, at soma level. The population LFP (p{t) is a 
sum of contributions 4>i{t) from cells at distances r^. The dependence of its amplitude a{R) on the 
population radius R serves to define the spatial reach (see text). B. Reconstructed cell morphologies 
used in simulations, L3 pyramidal cell (red), L4 stellate cell (green), L5 pyramidal cell (blue), 
superimposed on the layers' boundaries. Electrode positions shown as black dots. Brackets mark the 
regions where synapses were distributed for either homogenous (solid brackets), apical (dashed 
brackets), or basal (dotted brackets) stimulation. Adapted from [34]. 



are included, i.e., we study how the root mean square amplitude a of the population LFP (j){t) (obtained 
as a sum of single-cell contributions (j>i{t)) depends on the radius R of the subpopulation of cells included 
in the sum. 

In the simulations we use three different morphologically-detailed cell models shown in Figure [2|3: 
the layer-3 and layer-5 pyramidal cells, and the layer-4 stellate cells. All neuron models are passive, i.e., 
without active conductances, and the extracellular signatures of action potentials (spikes) are thus not 
included. For each class of pyramidal cells we consider three different spatial patterns of synaptic input: 
the synapses are placed either in the apical region only, in the basal region only, or evenly over the whole 
cell (Figure |2|3) . For the layer-4 stellate cells we consider only spatially homogeneous synaptic input, as 
these cells lack clearly defined dendritic regions. 

The synaptic currents are modeled as a- functions with a very short time constant (t=0.1 ms) to assure 
that no frequency filtering is imposed by the synapses themselves. In the frequency range considered in 
the present simulations (up to 500 Hz) each synaptic input current thus effectively corresponds to a S- 
function with a white (flat) power spectrum. With Poissonian spike statistics, which also implies a white 
power spectrum, the only frequency filtering in our simulation setup will come from the intrinsic dendritic 



filtering effect 40 42j due to electrical properties of the cable and the summation of the single-neuron 
LFP contributions to form the population LFP. For further details on the simulations we refer to the 
Methods section. 

Note that we here for simplicity will refer to all calculated extracellular potentials as "LFPs" even if 
we consider frequencies as high as 500 Hz, i.e., frequencies often regarded to be outside the LFP band. 

2.2 Simplified model of population LFP 

To understand how the population signal emerges from single-cell contributions we use a simplified 
mathematical model, which is a frequency-resolved version of the model introduced in [34]. 

We assume that the power spectral density (PSD) of the contribution to the LFP from the i-th cell 
at given frequency can be factorized as 

\Mf)\'^<^lif)Fnf), (2) 

where $j is the PSD of the single-cell LFP, cr| is the PSD of the input current, and Fi{f) is the frequency- 
dependent shape function of the i-th cell, which carries the information about how the root mean square 
amplitude of the signal at given frequency decays with distance. Moreover, we assume that the shape 
function of each cell in the population can be replaced with a single, distance- and frequency-dependent 
function: 

F,if) = F{f,r,), 
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that is, we assume that the shape function Fi only depends on the frequency and the lateral distance 
from the recording electrode (Figure [2|3), and neglect variation in the single-neuron LFP contributions 
due to other factors. For each particular morphology (layer-3/layer-4/layer-5) and synaptic stimulation 
pattern (homogeneous/apical/basal), the LFP contribution from each cell in the population is thus de- 
scribed with the function F{f, r). Note that for the special case of white-noise input (i.e., cr'^if) = const.), 
the squared shape function F{f,ri)'^ will be proportional to the PSD of the single-cell contribution to 
the LFP. 

The summation of single-cell LFPs to the population signal depends on the correlation between the 
single-cell LFP contributions. In the case of uncorrelated input this amounts to simply adding the 
variances of the single-cell LFPs. For a disc-like population of radius R we thus obtain the following 
expression for the PSD of the signal at the center: 



Go(/,i?)=a| J2 \F{f,n)\^ ^2alTTp T r |F(/, r)|' dr. 



(3) 



where p is the planar cell density. On the other hand, if the single-cell LFPs are fully correlated, the PSD 
of the signal is found by adding the single amplitudes, not variances, and we thus obtain 



ri<R 



rF{f, r)dr 



(4) 



In our simulation setup the single-cell LFP contributions from two equidistant neurons (i.e., same r^) 
are not identical even for Cin — 1: while the same spike trains are used to synaptically stimulate the cell, 
they will not in general activate an identical set of synapses (see Methods). Moreover, as we now work in 
the frequency domain, the correlation between single-cell contributions to the LFP (^i, 0j) is naturally 
replaced by their coherence ($*'i>j/|<i>i||<i>j|), which, in general, depends on the frequency. 

If we approximate the LFP coherence between each pair of cells by the population-averaged LFP 
coherence c$, then the PSD is given by 

P(/, R) = [l- Mf)]Go{f, R) + Mf)Gi{f, R) (5) 

(see Methods for the full derivation of this formula) . Note that the root mean square amplitude cr of the 
signal (see Figure [2]) is related to the PSD P{f,R) through 

a'{R)^ J P(/,i?)d/, 

where the integration is between / = Hz and half the sampling frequency. 

2.3 Numerical evaluation of ingredients of simplified model 

Equation [5] implies that any frequency dependence of the population LFP (for example, frequency depen- 
dence of the spatial reach) in general will result from the interplay of two separate effects: (1) frequency 
dependence of the single-cell shape functions F{f, r) and (2) frequency dependence of the coherence c,j,{f) 
between single-cell contributions to the population signal P{f, R). These two effects are addressed next. 

2.3.1 Frequency dependence of shape function 

The power of the extracellular potential from a single neuron decays when we move away from the cell, 
and the rate of the decay depends on the frequency of the signal. In Figure [3jV we have plotted squared 
shape functions F{f, r)^ at the soma level for three selected frequency bands for the case with the layer-5 
cell receiving basal synaptic stimulation. We observe that the high-frequency LFP component decays 
faster with distance than the low-frequency component. This leads to the low-pass filtered power spectra 
seen in Figure [3j3 and is consistent with our previous observations of low-pass filtering in dendritic cables, 
i.e., the intrinsic dendritic filtering effect |40[[42 . To quantify this phenomenon we approximate the actual 



shape functions with simplified power-law shape functions. Specifically, at the soma level the amplitude 
of the single-cell LFP is, following [45], modeled as: 

{Fo, if r < r^, 

Fo^/t^, ifr, <r<r,(/), (6) 

FoVrJr4f)ir4f )/r)^ if r > r,(/), 
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Figure 3. Ingredients of the simplified LFP model for the layer-5 cell with basal synaptic 
input. A. Spatial decay in lateral direction for the squared single-cell shape functions \F{f,r)\'^ for 
three different frequencies / = 0, 60 and 500 Hz. B. Single-cell LFP spectra \F{f,r)\'^ for three different 
lateral distances from the soma (dotted vertical lines in A). C. Log-log plot of the squared near-DC 

Hz) shape function |F(0,r)^| (dots) approximated by a piecewise-linear function with cutoff 
distance r* (line; see Eq. [6]). D. Frequency dependence of the cutoff distance ?■*(/). E. 
Population-averaged LFP coherence c$ for different input correlation levels Cin- Dots not connected 
with lines indicate that |c$| is plotted in place of spurious negative values (see Methods). F. Power 
spectra P(/, R) of the compound LFP (R = 1 mm); dots correspond to simulation; lines correspond to 
predictions from simplified model, Eq. [Sl based on r^, and c$ given in D and E, respectively. 



i.e., the shape function is approximated by oc r^^/^ close to the cell (r < r*) and by cx r^^ (dipole) in the 
far-field regime (r > r^). The parameter thus represents the cutoff distance where the LFP contribution 
switches from the near- field {F cx r-i/2) to the far-field regime (F (x r ^), see fitted curve in Figure pC. 
This parametric representation of the shape function allows us to express the functions Go(/, i?) and 
Gi{f,R) (Equations |3] and |4]) explicitly in terms of the cutoff distance r*, see Methods for details. The 
observed reduction of with increasing frequency (Figure [sjl)) is intimately related to the corresponding 
reduction of the frequency-dependent electrotonic length constant in dendrites jTO 40 . In the example 



shown in Figure [3|V the transition to dipole decay occurs closer to the cell for the high-frequency signal 
(at about ~ 40 /im) than for the low-frequency components (~ 80 /im). 

In Figure [4] we show the calculated cutoff distance r.^, for the various situations considered in the present 
paper involving the layer-3 pyramidal neuron (|4|\), the layer-5 pyramidal neuron (|4|3), and the layer-4 
stellate neuron (|4p). For the pyramidal neurons we consider three spatial patterns of synaptic inputs, 
that is homogeneous, only apical or only basal (Figure [2]). All these combinations of cell morphology and 
stimulation pattern exhibit similar behavior as in our example (Figure [3| : (/) decays with increasing 
frequency. The only exception is the layer-5 cell with apical input, where is very large, and also exhibits 
a minimum around 150 Hz. This reflects that the geometry of this situation is unique, with the synaptic 
input positioned far above the soma level where the LFP is recorded. As a consequence the shrinkage 
of the current dipole with increasing frequency will be accompanied by a shift of the mean position of 
the current dipole in the apical direction [34^ , 42^ . The squared shape functions and the single-cell power 
spectra for the remaining situations (all apart from layer-5 cell with basal synaptic input) are shown in 
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Figure 4. Frequency dependence of the cutoff distance for all situations considered: 

homogeneous (solid) , apical (dashed) and basal synaptic input (dotted) applied to the layer-3 pyramidal 
cell (A), the layer-5 pyramidal cell (B), and the layer-4 stellate cell (C). Cell morphologies depicted in 
Figure [2p. 



Figures S2A, B-S7A, B. 

2.3.2 Frequency dependence of coherence 

The single-cell shape functions F{f,r) alone are generally not sufficient to predict the population LFP. 
The missing component is cq,{J), the frequency-dependent population-averaged coherence between single- 
cell LFP contributions. This quantity can be estimated from population simulations, as described in 
detail in Methods, Equation |16[ Coherence curves for different input correlation levels for our example 
(layer-5 cells receiving basal stimulation) are shown in Figure jsjil. The coherence c$(/) is seen to be 
higher for low-frequency components. This may be understood on biophysical grounds by considering the 
dendritic morphology of the cell: for high-frequency synaptic input the return currents will be closer to 
the synaptic currents [40| , and for the example in Figure jsjil with basal stimulation of layer-5 pyramidal 
neurons, the resulting current dipoles will expectedly tend to be oriented randomly in space. However, 
for low-frequency input some of the synaptic input current will return through the apical dendrite f42j, 
and the orientation of the effective current dipoles will be more similar between cells, leading to a higher 
coherence. 

By combining the shape functions F{f,r) with the LFP coherence c^{f) in the simplified model 
(Equation [s]) we can now obtain predictions for the population LFP. The resulting PSD for our example 
situation is shown in Figure [sji^ and is seen to be in excellent agreement with the simulation results. (See 
Figures S2C-S7C for the results for the remaining combinations of cell type and synaptic input patterns.) 

In Figure |5] we show the frequency dependence of the coherence c$ (/) for the same full set of seven 
situations as depicted in Figure |4j A first observation is that for pyramidal neurons (layer-3, layer-5) with 
asymmetric synaptic input (either only apical or only basal), decay of c$(/) with increasing frequency is 
observed for all non-zero levels of input correlations Cin. This low-pass filtering effect is seen to be strongest 
for the layer-5 cell with basal input (Figure |5jV,[5|3,[5|l),[5ji;). However, when the same pyramidal neurons 
receive homogeneous synaptic inputs, the filtering effect is almost absent (Figure [5p,[5|i^). In that respect 
it resembles the situation with the stellate layer-4 cells receiving homogeneous synaptic input (Figure |5]GI) 
where c$ is essentially zero, implying that the correlations in the synaptic input do not translate into 
correlations of the single-neuron LFP contributions. 

2.4 Population LFP 

As a first step towards exploring the spatial reach of the extracellular potential in our disc-like setup we 
next show how the population signal emerges from single-cell contributions and investigate frequency- 
related effects. In Figure [6] we present results both from the full simulation and the simplified model 
(Equation |5| for our example situation with the population of layer-5 cells receiving basal synaptic input. 

In Figure |6}V we show the PSD of the LFP produced by differently-sized populations of cells receiving 
uncorrelated synaptic input. While we observe some low-pass filtering (especially above ~ 100 Hz) for 
all population sizes, the effect is not particularly strong. Figure [6|l) instead shows the PSD for the same 
uncorrelated situation as a function of the population radius R. We observe that the LFP in all frequency 
bands saturates rather quickly with increasing population size, that is for R ~ 100-200 fiui. This implies 
that the contributions from uncorrelated neuronal LFP sources positioned more than a few hundred 
micrometers away from the electrode are negligible for all frequencies considered. 
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Figure 5. Frequency dependence of the population-averaged LFP coherence c$ for all 
situations considered. Dots not connected with lines indicate that |c$| is plotted, see Methods. 



The situation changes dramatically for the case of correlated synaptic input (Figure [6|3,|6p,[6|5,[6|^), 
both in terms of amplitude and frequency dependence. For the case with the maximum input correlations 
Cin = 1 (Figure [6p, ^^), we see that the low-frequency power is up to two orders of magnitude larger 
than for the corresponding uncorrelated case. Further, a significant low-pass filtering effect is seen. For 
example, the low- frequency power (~ Hz) is an order of magnitude larger than the power at 60 Hz 
for Cin = 1 (Figure [gJ^"). Another observation is that the low- frequency power grows much faster with 
increasing population radius than the high-frequency power (Figure [6|5, [6j^). Finally, the power of the 



population signal no longer seems to saturate as the population radius increases 34 

The predictions from the simplified model agree qualitatively with the full simulation results; however, 
we observe some clear deviations: First, in Figure [6)l)-F we see that the simplified model overestimates 
the power of the low- frequency components (~ Hz, 60 Hz). This is because the model here uses the 
approximate power-law shape functions (Equation [6| which lie above the numerically evaluated shape 
functions for low frequencies (Figure [3p). For high-frequency components (500 Hz), on the other hand, 
the opposite situation occurs (results for fitted approximate power-law function not shown). Second, in 
case of correlated input the model works better for the larger populations than for smaller ones. This is 
as expected given the present procedure for calculating the LFP coherence c$(/) used in the simplified 
model: here this LFP coherence c$(/) was extracted from the full population {R = 1000 /im) simulations, 
and the value obtained is not surprisingly a poor approximation when applied to populations which are 
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Figure 6. Power spectral density of population LFP as a function of frequency and the 
population radius. Full simulation results (dots) and simplified model predictions (lines) for the LFP 
at the center of disc- like populations of layer-5 pyramidal cells receiving basal synaptic input. Three 
different input correlation levels cin are considered. A-C: PSD of population LFP for three population 
radii R. D-F: dependence of power of three different frequency components on the population radius R. 



much smaller. With c$(/) calculated for each population radius R separately, the simplified model 
predictions significantly improve (Figure SI). 

2.5 Frequency-dependence of spatial reach 

We are now ready to analyze the frequency dependence of the spatial reach of extracellular potential. 
Following |34| we define the spatial reach as the radius of the subpopulation which yields 95% of the 
root mean square amplitude in the population center compared to the largest population considered 
{R = 1 mm). With this definition the spatial reach is easily found from the data presented in Figure [6p, 
[6^5 andlGj^" as the distance at which the amplitude of the LFP reaches 95% of the maximum value. 

The results for the spatial reach for all seven situations considered are shown in Figure [7] The reach 
is seen to vary both with the frequency / and the level of input correlation Cin, but the specific effects 
depend sensitively on the cell morphology and synaptic stimulation pattern. For the pyramidal cells with 
asymmetric input (either only basal or only apical) the spatial reach grows significantly with increasing 
input correlations Cjn (Figure [7|V,[7|3.[7p,|7|i;). The effect is particularly prominent for lower frequencies, 
i.e., smaller levels of input correlations Cin are needed to increase the spatial reach significantly. As 
a consequence, for certain correlation levels Cjn the spatial reach of the low-frequency components can 
differ a lot from the spatial reach of the high-frequency components. For example, in the situation 
with the layer-5 population receiving basal input with Cjn = 0.01, the spatial reach at 100 Hz is only 
around 200 fim, while the low-frequency reach is almost 800 /im. For the case of homogeneous inputs 
into pyramidal neurons (Figure |7p, [tJ^") these effects are still present, but seen to be much weaker. For 
the layer-4 stellate cells the spatial reach is practically independent of the frequency / and the input 
correlation level Cjn, Figure |7pl. 

Note that the situation with the layer-5 population receiving only apical input is again somewhat 
different from the other cases. Here the spatial reach for the uncorrelated input is already quite large 
(~ 300-400 /im) and the levels of the input correlation required to saturate the spatial reach at a maximum 
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Figure 7. Spatial reach for different frequency components of LFP. Spatial reach is defined as 
the radius of a subpopulation contributing 95% of the root mean square amplitude of LFP compared to 
the amplitude for R — 1000 /im. LFP is calculated at the soma level at the center of the population. 
Full simulation results plotted with dots; predictions from the simplified model (Equation [s]) based on 
calculated values of r^, and c$ given in Figures l4| and [Sl respectively, are shown with lines. 



value possible in our setup are significantly smaller. 

For the case of uncorrelated input we can obtain analytical expression for the spatial reach from the 
simplified model. Using Equations [s] and [6] we obtain an explicit formula for Go{f,R) in terms of the 
cutoff distance r^,(f), the population radius R and r^. From this, we find in the limit of — ^ 0, that the 
radius of the subpopulation contributing a fraction a of the asymptotic amplitude (i? — > cxd) is equal to 
r^/VS — 3q;2 (valid for > |) . For our choice of a = 0.95 we find the spatial reach to be ~ 1.85r*. 

2.6 Decay of extracellular potential outside the population 

The spatial reach we have discussed above represents an 'electrode-centric' point of view: we ask about 
the distance from the recording electrode of the neurons setting up the LFP signal. However, one can also 
take a 'population-centric' approach and instead ask how rapidly the LFP signal decays with distance 
outside an active population [34'. 

In Figure [8] we show results for this situation for an example population of layer-5 cells receiving 
basal or apical synaptic inputs. The first observation in case of basal synaptic input is that the low- 
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and medium-frequency LFP components Hz, 60 Hz) are significantly boosted, up to two orders of 



8 3) . This applies both inside and outside of 
), however, input correlations are seen to 



magnitude, by high levels of input correlations Cjn (Figure[8j/ 
the population. For the high-frequency signal (500 Hz, Figure 
have only a small boosting effect on the signal amplitude. In the case of apical synaptic inputs the effect of 
increasing input correlations is seen to be more uniform across frequency bands, with the high-frequency 
components (500 Hz) being boosted by roughly the same factor as the low- and medium-frequency LFP 
components Hz, 60 Hz), Figure P 




Figure 8. Decay of extracellular potential outside populations of layer-5 cells with 
asymmetric input. Each of the panels shows full simulation results (dots) and predictions from 
simplified model. Equation [S] (lines) for one frequency band (0, 60, 500 Hz) and four input correlation 
levels. Horizontal dotted lines indicate 'noise level' (power of the signal generated by a population of 
uncorrelated cells with homogeneous input, see text). A, B, C: basal synaptic input. D, E, F: apical 
synaptic input. 

The strong boosting of the LFP signal seen for correlated synaptic input for ^ Hz (Figure [8|V) 
and 60 Hz (Figure [8|3) has direct implications for how recorded LFP signals should be interpreted. 
As observed in these panels, highly-correlated populations some distance away from the electrode may 
easily dominate contributions from uncorrelated populations surrounding the electrode. For example, 
in Figure |8jV we observe that the LFP signal 500 /im outside a correlated population with Ci,i=0.1 is 
still larger than the contribution recorded inside the same population receiving uncorrelated synaptic 
inputs (cin=0). For 60 Hz (Figure [8|3) the boosting effect is smaller, but still the signal recorded outside 
a correlated population may be larger than what is recorded inside an identical population receiving 
uncorrelated input. The same effect is seen to be even more pronounced for the apical-input case in the 
lower panels (Figure [8|l)-[8|^), further highlighting that the interpretation of the recorded LFPs in terms 
of activity in the neurons immediately surrounding the electrode has to be done with caution. 

In Figure [9] we show the same PSDs as in Figure [8] but normalized to unity at the population center. 
This illustrates that the decay of the LFP is more abrupt around the population edge in the uncorrelated 
case than in correlated cases (this is especially prominent for the low-frequency components ^ Hz, 
60 Hz). This is consistent with an observation made in |45| (see Figure 3.9 therein), namely that in the 
large-population limit the LFP signal power at the population edge will be reduced to half of power at 
the center for uncorrelated populations, while it will be reduced to a quarter of the center power for 
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fully correlated populations. Here this difference between the correlated and uncorrelated cases is more 
pronounced for the low-frequency components, where the coherence c$ is largest. 
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Figure 9. Decay of extracellular potential outside populations of layer-5 cells with 
asymmetric input. Same as Figure [8] but with PSDs normalized to 1 at the population center, and 
the distance axis zoomed in to highlight the behavior around the edge of the population. 

We next investigated the related question of detectability, i.e., how far away from a synaptically 
activated population the generated LFP still can be detected above the ambient LFP 'noise'. This noise 
level will naturally vary between experimental situations, but here we assumed it to be given by the 
background LFP signal from neurons receiving the same number and type of synaptic inputs, except 
that (1) the inputs are uncorrelated and (2) homogeneously spread over the neuronal membrane. (The 
power of this background LFP signal is plotted as dotted lines in Figure |8]) The frequency-dependent 
signal decay and detectability outside basally-activated populations are illustrated in the 2D color plots 
in Figure [T0| As in Figure [8] the population radius is fixed at i?=1000 ^m, and we plot the PSD both 
inside and outside the population. The lines mark where the signal-to-noise ratio falls below 0.5 (solid 
line) and 0.1 (dotted line), respectively. Here the signal-to-noise ratio is defined as the ratio between 
the root mean square amplitudes of the LFP signal (from the basally-activated population) and the LFP 
noise (from the background population). 

A first observation is that for uncorrelated synaptic inputs (cin=0. Figure 10 \.-[To)3), there is very 



little variation with frequency. Also the detectability of the LFP outside the active population is poor: 
the signal-to-noise ratio falls to 0.5 about 100 /im outside the population, and below 0.1 less than 500 /im 
outside. The situation is seen to be very different when the populations receive correlated synaptic inputs. 



Focusing first on the case with the largest level of input correlations (cin=l. Figure 10 GI, 10 3), we see 
that the lower frequencies of LFP extend further outside the population than the higher frequencies. 
For example, for the near-DC component (~ Hz) the signal-to-noise ratio is seen to be almost 0.5 at a 
distance of 2000 /im, i.e., 1000 /im outside the population edge, and 0.1 as far way as 2000 /im outside this 
edge. For the 125-Hz component, on the other hand, the signal-to-noise ratio is reduced to 0.5 as little 
as 200 /im outside the population. The results for the intermediate cases (ci„=0.01, Cin=0.1) depicted in 
Figures [TOp-[lO)^ are seen to bridge these uncorrelated and strongly correlated cases. 

The results for the basally-driven pyramidal cell population in Figure [TO] demonstrate a main result 
from this study, namely that correlations in synaptic inputs may significantly enhance the amplitude 
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and thus also the detectabihty of the low-frequency LFP components relative to the high-frequency LFP 
components. The same effect is observed for the same population when the synaptic inputs are placed 
solely on the apical part of the neurons, cf. Figure [TT] However, here a sizable low-pass filtering effect in 
detectability is observed also for the case with uncorrelated input (Figure [TTJV,[TTJ3) due to the intrinsic 
dendritic filtering effect [40' 42' . It is also worth noting that populations of layer-5 cells stimulated apically 
yielded the farthest-reaching LFP signal of all cases analyzed. Note also that the low-pass filtering effect 
in the boosting of LFP signal with increasing correlations was seen to be largely absent in the case of a 
spatially homogeneous distributions of synaptic inputs onto populations made of any of our three example 
neuronal morphologies (results not shown). 

Finally, inspection of Figure [8]( and the PSD line plots in Figures 10 and 11 1 reveals that the predictions 
from the simplified model (Equation |5| agree excellently with the full numerical simulations for the case 
of uncorrelated input. However, the simplified model systematically overestimates the signal power for 
correlated populations for positions far outside the active populations. This is because the simplified 
model predicts a fall-off of the LFP amplitude proportional to in the far-field limit, while in the 
full simulations the total LFP signal will be dominated by correlated dipoles oriented vertically. As a 
consequence the functional form of the lateral signal decay will be closer to (34]. 

This limitation of the simplified model can be remedied by incorporating the fact that the evaluated 
population-averaged coherence C(j>(/) not only depends on the size of the population R considered, but also 
on the electrode position X along the horizontal axis from where it is evaluated, i.e., c,j,{f) = Cii,{f; R, X). 
So far the population-averaged LFP coherence has been evaluated at the population center, i.e., at X = 0, 
but when Equation [16] is evaluated at other positions X, as shown in Figur e |12[ c$ is observed to decay 
as 1/ X^ for X ^ R. In the formula for the simplified model in Equation l5] the power P is seen to be 
proportional to c$Gi in the correlation-dominated regime. A modified simplified theory including not 
only the X-dependence of Gi [34j|45j, but also the observed X-dependence of c$, indeed predicts the 
correct far- field X-dependence outside the active population (see Figure S8). The physical interpretation 
is that the dominance of the LFP signal of the correlated vertical dipoles will be incorporated in the 
population-averaged LFP coherence c$. 



3 Discussion 

In this computational study we have investigated the frequency dependence of the signal power and 
'locality' of cortical local field potentials (LFP). While some low-pass filtering effects of the LFP are 
seen also for populations of cells receiving uncorrelated synaptic inputs or homogeneously distributed 
correlated synaptic inputs, the large frequency-dependent effects are seen when populations of pyramidal 
neurons receive correlated and spatially asymmetric inputs (i.e., either only basal or apical). For example, 
for the case with a layer-5 population receiving correlated, Poissonian synaptic currents (with a white- 
noise, i.e., flat band, power spectra) onto their basal dendrites, the power of the low-frequency LFP 
(~ Hz) was seen to be an order of magnitude larger than the corresponding LFP power at 60 Hz. 
Correspondingly, the low-frequency LFP components was seen to extend much further outside the active 
population than high-frequency components. 

The findings from our comprehensive biophysical simulations using reconstructed neuronal morpholo- 
gies were backed up by a simplified model, adapted from |34 , for generation of population LFP. This 
model is based on three factors: (1) the decay of the single-neuron contribution with the distance from the 
electrode represented by the frequency-dependent shape function F(f^r), (2) the population geometry 
and density of neuronal LFP sources, and (3) the frequency-dependent correlation (or, more precisely 
coherence c$(/)) of the single-neuron LFP contributions from individual neural sources. Our simple 
model for the population LFP (Equation |5| was found to give quantitatively accurate predictions, im- 
plying that it captures the salient features. While some of the observed low-pass filtering could be traced 



back to single-neuron properties and the intrinsic dendritic filtering effect 40 42' accounted for by the 
shape function F{f,r), most of the observed low-pass filtering was due to strong low-pass filtering in 
the coherence c^{f) between the single-neuron LFP contributions: synaptic-input correlations translated 
into correlated single-neuron LFP contribution to a much larger extent for lower frequencies than for 
higher frequencies. As a direct consequence, the low-frequency components of the extracellular potential 
are significantly boosted in populations with correlated synaptic input. In our model this happens purely 
because of dendritic filtering, as the synaptic input currents themselves have been tailored to have a 
fiat (white-noise) PSD. With a colored (frequency-dependent) spectrum of the synaptic input, the power 
spectrum of the LFP would be given as the product of the PSD of this synaptic filter and the PSD from 
the dendritic filtering investigated here. 

A key qualitative finding in our study is that the size of the signal-generating region, i.e., the spatial 
reach, may in the case of correlated synaptic input vary strongly with frequency. For the example 
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Figure 10. LFP signal power as functions of frequency and distance from basally-activated 
pyramidal-cell populations. Colormaps (A, C, E, G) show the power of extracellular signal of a 
population of layer-5 cells receiving basal synaptic input for four levels of input correlation Cin as 
functions of frequency and distance from center of populations. Black solid and dotted lines denote 
signal to noise ratio of 0.5 and 0.1, respectively. B, D, F, H: power spectra of extracellular signal at 
different distances, lines: prediction from simplified model in Equation |5] dots: full simulation. Thin 
vertical dotted lines with dots in A, C, E, G denote the distances at which the power spectra are shown, 
that is, at the center (0 /im), population edge (1000 fim), and two distances outside 1600 /im and 
- 2500 fnu). 



population in Figure [TJwc sec that for Ci„ = 0.01, a plausible correlation level in cortical spiking networks 
(see, e.g.. Figure 6 in |34|), the LFP spatial reach may be reduced from close to the size of the population 
(^ 800 /im) for ^ Hz to ^ 400 /im for 60 Hz. For uncorrelated input, however, the spatial reach 
will generally always be small (< 200 /im) for all frequencies, with the exception of the case with apical 
input on large pyramidal cells (Figure [7| . Note that in the present simulation scheme the spatial reach is 
by definition less than 1000 /im, the size of our model population. Unlike for uncorrelated populations, 
the LFP power will for correlated populations keep on increasing when the population grows beyond 
1000 /im |34j. The present definition of spatial reach (95 % of the amplitude for i?=1000 /im) thus 
underestimates the true size of the signal-generating region in this case. 

In a recent experimental study from macaque auditory cortex [33' it was observed that different 
frequency bands spread equally far from a source (cf. Figures 5 and 6 there). There are, however, notable 
differences between this study and our present approach, making it difficult to compare the results. First, 
here we investigate the spread of the LFP along cortical layers at the soma level, while in fSS] the spread in 
vertical direction was studied. Second, and likely more importantly, in [33^ the LFP amplitude at a given 
latency after stimulation was used to extract LFP decay profiles. In contrast, we here use noise input 
and consider the root mean square amplitude of LFP over a relatively long time period. Further, the 
correlation level of the synaptic input, found here to be a critical parameter in determining the frequency 
dependence, is not known in the situation in |33 . It is thus difficult to assess whether our results are in 
accordance, or not. 

Our results have direct consequences for the interpretation of observed cross-correlations between 
extracellular potentials recorded at different electrodes |46-51 . As demonstrated here the low-frequency 
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Figure 11. LFP signal power as functions of frequency and distance from 
apically-activated pyramidal- cell populations. Same as Figure 10 but for a population of layer- 
cells with apical synaptic input. 



LFP signal generated by a population of neurons around one electrode receiving asymmetric synaptic 
input, may extend a millimeter or more outside the active population (see, e.g.. Figure 10 GI). Thus 



measured correlations in the low- frequency LFP components between two electrodes positioned, say, one 
millimeter apart, may be due to volume conduction effects. However, cross-correlation induced by such 
volume conduction will, as demonstrated here, have a diminishing spatial range with increasing LFP 
frequencies. Note also that the magnitudes of the LFP amplitude at the two adjacent electrodes will 
aid in the interpretation: while volume conduction may propagate the LFP a millimeter or more, the 
amplitude will rapidly diminish with distance (cf. Figure 10 and 111. Thus the observation of large- 
amplitude LFPs at both electrodes would be an indication that both electrodes are surrounded by strong 
LFP-generating populations. 

In [52 the temporal power spectra of the EEG were shown to be well fitted by 1//" power-law 
functions with power-law exponents a varying between brain areas: in the frontal lobe a was reported to 
be 1.78 ± 0.76, while in the occipital lobe a — 1.19 ± 0.28. Power laws have also been found in recordings 
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Figure 12. Population-averaged LFP coherence c$ as a function of distance X from center 
of population of layer-5 pyramidal cells. A, B, C: basal synaptic input, D, E, F: apical synaptic 
input. Dots not connected with lines indicate that \c,^\ is plotted in place of spurious negative values 
(see Methods). Dashed lines mark decay. 



of the LFP, see, e.g., 53 54 , often with different exponents a. In [54] a was shown to vary between 
network states, more specifically between the slow-wave sleep and awake states. In this context it is 
interesting to note that the PSDs in our Figure express approximate power laws with exponents a 
highly dependent on the degree of coherence. This finding suggests that varying levels of coherence in the 
synaptic input may be a mechanism underlying the different experimentally observed power laws. This 
would also be in agreement with the experimental observations that network states with a presumably 
large coherence (e.g., slow wave sleep in [54]) typically express a larger value of a than network states for 
which the coherence is lower (e.g., awake state in [54]). 

In our modeling we have assumed the extracellular medium to have a frequency-independent con- 
ductivity, an assumption supported by a recent thorough experimental study of the electrical properties 
of monkey cortical tissue [55^ . However, if for example low- frequency filtering (Tcond = '^mndif) of the 
extracellular medium should be found [56] , this filtering would superimpose directly on the filtering seen 
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here, i.e., the total LFP filter would be the product of the LFP filter calculated here and the filter from 
the extracellular medium (oc l/<Tcond(/))- 

Here we have focused on the spatial and spectral properties of LFP signals triggered by presynaptic 
spikes that could originate from within the same cortical population or come from other distant brain 
regions. While not addressed here, it may be that the LFP signal itself influences the timing of these 



locally generated spikes through ephaptic coupling 57 58 . That would in turn influence the correlation 



structure of incoming spikes and thereby also the generated LFP signal. Since our simulations show that 
both the LFP amplitude and spatial reach is larger for low than for high frequencies, this suggests that 
if ephaptic effects play a role in cortical processing, they would likely be larger for low than for high 
frequencies. 

In the present analysis we have modeled the dendrites as simple RC-circuits which, in combination 
with the use of current- based synapses, made the system linear. This greatly facilitated the present 
frequency-resolved analysis in that the LFPs at different frequencies were effectively decoupled, cf. the 
standard theory for Fourier analysis of lineax systems. The present results also serve as a starting point 
for the exploration of non- linear effects, for example due to active membrane conductances. Close to 
the resting potential of the neuron, the active conductances can be linearized, and the neuron dynamics 
can be described by linear theory with quasi-active membrane modeled by a combination of resistors, 
capacitors and inductors (see, e.g., Ch. 10 in [59], Ch. 9 in [60], or (6l]). At present it is not known to 
what extent such 'generalized' linear schemes will be able to account the LFP generation in real neurons, 
but the present forward-modeling scheme, applicable for passive and active conductances alike, can be 
used to explore this question systematically. 



4 Methods 

4.1 LFP simulations 

The setup of the LFP simulations is almost identical to the scheme used to model cortical population 
LFPs in ^34] . The main difference is that here we use a much smaller synaptic time constant to achieve 
an effectively white (flat) power spectrum for the synaptic currents for the frequencies of interest here 
(less than 500 Hz). We therefore also use a smaller numerical time step. The model parameters are 
presented in detail (in the format described in |62]) in Tables SI, S2 and S3. For the reader's convenience 
we summarize the essential information below. 



4.1.1 Cell models 

We analyze three conipartmental cell models: the layer-3 and layer-5 pyramidal cells, and layer-4 stellate 
cells [63], available from ModelDB fO?, accession number 2488. We modified the models by removing 
active conductances and axon segments. The passive parameters of the cells were the following: spe- 
cific axial resistance Ra — 150n-cm, specific membrane resistance R„i = 30 kfi-cm, specific membrane 
capacitance Cm = 1.0/iF/cm. 

Each simulated cell was stimulated using 1000 excitatory current-based a-function synapses with a 
time constant t = 0.1 ms. The synaptic time constant was short enough to ensure that the spectrum 
of the input current was flat in the studied range. Each synapse was driven by a homogeneous Poisson 
spike train with the rate of 5 spikes per second. The spike trains driving one cell were independent. For 
uncorrelated input into the population also the spike trains belonging to each cell were independent, for 
correlated input they were drawn (without repetitions for each cell) from a common pool consisting of 
1000/cin spike trains. As a result, in case of correlated input each two cells shared 1000 • Cin spike trains 
on average. Note that even for Cjn = 1, when each of the cells is driven by the same spike trains, the 
spike trains will in general be assigned to different synaptic locations. 

We simulated activity of cells for either 10200 ms (single-cell shape functions and LFP in the pop- 
ulation's center) or 1200 ms (LFP at points not in the population's center). The first 200 ms were 
discarded to avoid start-up artifacts. We used a fixed time step of 1/64 ms, and recorded the results of 
the simulation (transmembrane currents in all compartments) with 1 ms time step (sampling frequency 
1 kHz). 

For the pyramidal cells we employed three stimulation patterns: the synapses were distributed either 
in the apical or basal part, or homogeneously throughout the whole dendritic tree (in each case the 
probability of attaching a synapse in a given compartment was proportional to its surface area). We used 
the same layer boundaries and soma depths as in [34]. 
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4.1.2 Calculation of LFP 

The extracellular electric potential was calculated using the line-source method [37 65 , resulting from 
integration of Equation [l] over linear dendritic segments. We assumed a purely resistive, homogeneous, 
isotropic and infinite extracellular medium, and an ideal point electrode (no filtering) , placed at the soma 
level. In single-cell simulations the electrode was placed at a distance (between 10 /im and 10000 fim) 
from a single cell, in population simulations it was placed either at the center of the population or at 31 
points placed between fim and 10000 iim from the center. To obtain the model LFP at the center of 
differently-sized populations we summed contributions from different subsets (cells located closer to the 
electrode than some distance) of the same full {R = 1000 /im) population. 



4.1.3 Single-cell shape functions 

To obtain single-cell shape functions (Figure |3jV) we calculated the LFP at different distances from a 
single cell, then calculated power spectra of these signals. The final curves were obtained by averaging 
power spectra from 100 simulations for each distance. 



4.1.4 Population simulations 

We simulated populations of = 10000 identical neuron. The cells were placed homogeneously within a 
disc of 1 mm radius at the same depth. Each cell was rotated randomly along the vertical axis. 



4.1.5 Software 

We performed the simulations using the NEURON simulator ( |44| , www.neuron.yale.edu) and the 
Python (www.python.org) interface to NEURON f66'; we also used NeuroTools (neuralensemble.org/ 
trac/NeuroTools/). The calculations of extracellular field were performed using LFPy |67| — Python 
package for modeling of LFP. 



4.2 Derivation of the mean-field model 



To derive the formula in Equation [s] for the power spectral density (PSD) of the extracellular signal in 
the center of the population we start we the assumption that the PSD of the contribution of the i-th cell 
at given frequency / (<i>i(/)) may be factorized as 



mf)\'^^lif)Fi{f), 



(7) 



where is the PSD of the input current, and Fi{f) is the frequency-dependent shape function of the 
i-th cell. We also assume that the shape function F depends only on frequency and distance from the 
center, that is: 

F,{f)^F{f,n). (8) 

Let us compute the PSD of the population signal $(/) (dependence on frequency / dropped below 
for convenience): 
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N N 



p = i$(/)p = <f = E E 'J'^ = E 'f^**' + E E*^**. 



(9) 



1=1 j=i 



We now use Equations [7] and [8] to express P in terms of shape functions and the PSD of the input current, 
note the trick in the double sum: 



P = ai 



V 



(10) 



We further assume that the coherence term |^|^ may by replaced by its population average over 
N{N — 1) pairs, hence we can move it in front of the double sum: 



/ N N N 



^ ' i=l j = l I 'II J I i=l 7 = 1 



(11) 
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As marked in Equation [TT] we denote the population-averaged coherence by c$. We further rewrite P as 



N N ^ N 



i=l 



i=l 



and finally 



N 



N 



P = (1 - c$) a| ^(^«)' +c* 4\ ^(^») 



(12) 



(13) 



^ — 



If we assume approximate, power-law shape functions F(/, r) parametrized by the cutoff distance 
f{J) (Equation [6]) , and change sums to integrals as in Equations |3] and [4] (limit of large number of cells), 



then the functions Go{f,R) and Gi(/, i?) have the following closed-form representation 45 



Go{f,R) 



Gi(/,i?) = 



a R<r„ 
if < i? < r* , 
if r, < R, 



F^PttR^ 

F§pTrr,{2R~r,) 
F^p^r,{ir,-r,-rl/R^ 

'F^P^tt^R^ 
lF§p^7r^r,(r!/^ - {4 + 6\og{R/n))r'J^ 



ifR<r„ 

if re < i? < r* , 

if r, < R, 



(14) 



(15) 



which we used for calculating predictions from the simplified model. The model can be modified to 
calculate the power of the signal outside the center of the population, i.e., at positions offset from the 
center by the distance X. In that case, the function F{f, R) in ^ and ^ has to be replaced by 
F{fj k ~ ^D- It is no longer easy to obtain closed-form formulae for Go and Gi in terms of r*, and we 
used the (non-parametric) shape curves obtained from the simulations, as the final integration had to be 
done numerically anyway. 



4.3 Data analysis 

4.3.1 Population-averaged LFP coherence 

It is hard to estimate the population-averaged LFP coherence c<i> directly as an average of pairwise 
coherences between the single-cell contributions to the LFP. Therefore, we used the same technique as 
in (34] (Equations 14 and 15 therein), ending up with 



c*(/) 



iV(iV- 



N 



1) 



(16) 



Coherence is always positive; however, the population-averaged coherence c<j estimated using Equation |16| 
may take spurious negative value (for example because finite-length signals are used) . This does not mean 
that c$ is truly negative, but rather that the value is too small to be estimated reliably from the amount 
of data available. In such cases we plotted |c$| in figures. 

Note that *&i(/) in Equation 16 may be evaluated either at the population center, or at a lateral 
position X > 0; as a result we will get either c$(/; X = 0) or cg,{f; X), see the last Section of Results. 

4.3.2 Frequency analysis 

To calculate the power spectral densities we used the Welch's average periodogram method (the matplotlib . 



mlab.psd function from Matplotlib 68 ). We used a Hanning window of length 32 or 128 time steps (32 
or 128 ms) and overlap between blocks equal to the half of the window length, which resulted in 17 (or 
65) equally spaced frequency bins between and 500 Hz. When calculating the population-averaged LFP 
coherence, Equation [16] we calculated the discrete Fourier transform and binned the resulting c$ in the 
same frequency bins as resulting from the Welch's average periodogram method. 



4.3.3 Spatial reach of LFP 

To obtain the spatial reach of the LFP we calculated the power spectral density P(/, R) of the population 
LFP as a function of frequency / and population radius R (taking values between and 1000 /im in 25 fiiJi 
increments). The spatial reach at given frequency was defined as the smallest radius R* for which the 
amplitude (T,s>{f,R*) is larger than 95% of the amplitude calculated for the full population. 
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4.3.4 Software 

Data analysis was performed using NumPy and SciPy Python packages [69] and IPython \70j. Plotting 
was done using Matplotlib [68] , 
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Supplementary Information 

Frequency dependence of signal power 
and spatial reach of the local field potential 

Szymon L^ski, Henrik Linden, Tom TetzlafT, Klas H. Pettersen, Gaute T. Einevoll 





Figure SI. Power spectral density of population LFP as a function of frequency and the population 
radius. Full simulation results (dots) and simplified model predictions (lines) for the LFP at the center 
of disc-like populations of layer-5 pyramidal cells receiving basal synaptic input. Three different input 
correlation levels Cin are considered. A-C: PSD of population LFP for three population radii R. D-F: 
dependence of power of three different frequency components on the population radius R. This is an 
alternate version of Figure 6 from the paper; here the coherence c$ is estimated not just once for the 
full (R = 1000 ^m) population, but in a radius-dependent fashion, for each population radius 
R — 25, 50, 75, . . . 1000/im separately. In effect the simplified model predictions are closer to the full 
simulations than in Figure 6. 
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Figure S2. The shape function \F{f,r)\'^ and the population LFP power spectra for layer-5 cells with 
apical input. A. Spatial decay in lateral direction for the squared single-cell shape functions \F{f,r)\'^ 
for three different frequencies f = 0, 60 and 500 Hz. B. Single-cell LFP spectra \F{f,r)\'^ for three 
different lateral distances from the soma (dotted vertical lines in A). C. Power spectra P{f, R) of the 
compound LFP (R = 1 mm); dots correspond to simulation; lines correspond to predictions from the 
simplified model. 
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Figure S3. The shape function \F{f,r)\'^ and the population LFP power spectra for layer-5 cells with 
homogeneous input. See Caption of Figure S2. 
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Figure S4. The shape function \F{f,r)\'^ and the population LFP power spectra for layer-3 cells with 
apical input. See Caption of Figure S2. 
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Figure S5. The shape function \F{f,r)\'^ and the population LFP power spectra for layer-3 cells 
basal input. See Caption of Figure S2. 
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Figure S6. The shape function r)p and the population LFP power spectra for layer-3 cells with 
homogeneous input. See Caption of Figure S2. 
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Figure S7. The shape function \F{f,r)\'^ and the population LFP power spectra for layer-4 cells with 
homogeneous input. See Caption of Figure S2. 




Figure S8. Decay of extracellular potential outside populations of layer-5 cells with asymmetric input. 
Each of the panels shows full simulation results (dots) and predictions from simplified model Equation 5 
(lines) for one frequency band (0, 60, 500 Hz) and four input correlation levels. Horizontal dotted lines 
indicate 'noise level' (power of the signal generated by a population of uncorrelated cells with 
homogeneous input, see text). A, B, C: basal synaptic input. D, E, F: apical synaptic input. This is an 
alternate version of Figure 8 from the paper, here the population-averaged coherence c$ depends also on 
the distance. In effect the simplified model predictions are closer to the full simulations than in Figure 8. 
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A 


Model summary 


Population 


cylindrical homogeneous cortical populations 


Neuron 


passive multi-compartment neuron models 


Synapse 


current based, alpha-shaped postsynaptic current with short time constant 


Input 


uncorrelated/correlated Poisson spike train input or spike train input gener- 
ated by laminar network model 


Measurements 


simulated LFP 


B 


Population 


Type 


separate homogeneous populations consisting of N neurons 

Population types: L3 pyramidal cell population, L4 stellate cell population, L5 
pyramidal cell population 


Geometry 


cylinder of radius R ('cortical column') subdivided into layers (see Table S3A) 


Cell positions 


- random soma positions on a disc at soma depth in vertical midpoint of 
corresponding cell type k, soma density p = N /(ttR^) 

- random cell rotations along vertical cylindrical axis 


Parameters 


N , R, Zfe, layer boundaries 


C 


Neuron 


Type 


passive multi-compartment neuron models with reconstructed morphologies 


Morphology 


- L3 pyramidal cell 

- L4 stellate cell 

- L5 pyramidal cell 

from [63], downloaded from ModelDB, accession number 2488 
axon compartments were removed 


Neuron 
dynamics 


non-spiking neurons with passive membrane with specific membrane resistance 
i?,n, specific axial resistance R^, and specific membrane capacitance Cm 


Compartments 


segments length shorter than one tenth of electrotonic length for 100 Hz re- 
sulting in 549 compartments for the L3 cell, 343 compartments for the L4 cell 
and 1072 compartments for the L5 cell (for chosen passive parameters and 
morphologies, see Table S3) 


Parameters 





Table SI. Summary of population model used in LFP simulations. Continues in Table S2 
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D 


Synapse 


Type 


current synapse with a-function shaped postsynaptic currents (PSCs) 


Dynamics 


input current: I{t) = ^ - h{t — ti) 

with kernel h{t) = /o;f ei-(*-*')/^6l(t - U) 

where ti is the arrival time of presynaptic spike i, 9 is Heaviside step function, 
with time constant r and amplitude /o= I^, /j for excitatory and inhibitory 
synapses respectively 


Parameters 


^1 



Input 



Poissonian input 

Uncorrelated each cell receives risyn independent excitatory inputs generated by Poissonian 

input point processes with rate rs„„ 



Correlated 
input 



each cell receives risyn excitatory inputs drawn (without resampling) from a 
finite-sized pool (size ripooi) of independent Poissonian point processes with 
rate rgyn resulting in correlation Cin = 'T-syn/?^pooi between total input current 
of different cells 



Synapse 
placement 



synapses placed on dendrites in certain layers depending on synaptic input 
region (apical/homogeneous/basal) and cell type: 



cell type 


apical 


homogeneous 


basal 


L3 


upper half of L23 


LI and L23 


lower half of L23 


L4 




L4 




L5 


LI and L23 


all layers 


L5 and L6 



Random placement of synapses within allowed boundaries with uniform den- 



sity with respect to membrane area (note: 
region) 



is fixed irrespective of input 



Parameters 



''syn 3 ' syn i 



F Measurements 


Simulated LFP 


Type 


extracellular field potentials (representing the LFP) calculated using the line- 
source method 


Assumptions 


extracellular medium as assumed to be purely resistive (non-capacitive) and 
inifitely-volumed with extracellular conductivity (Tcond 


Electrode 
properties 


ideal (non-filtering) point electrode placed either in the center (r=0) of the 
cylindrical geometry, or offset by some distance; at depths corresponding to 
depth of somata (see Table Sl:Population). 


Parameters 


'^cond 



Table S2. Summary of the population model used for LFP simulations. Continued from Table SI 
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Name 


Description 


Value 


N 


number of cells 


10000 


R 


population radius 


1000 ^ni 



Zk and layer 
boundaries 



layer boundaries and soma positions (in units of /xm in relation to cortical surface 
at /um), derived from Stepanyants et al., Cereb Cortex 2008, 18(l):13-28: 



layer 


upper boundary 


soma depth Zk 


lower boundary 


LI 







-81.6 


L2/3 


-81.6 


-334.3 


-587.1 


L4 


-587.1 


-754.6 


-922.2 


L5 


-922.2 


-1021.1 (1) 


-1170.0 


L6 


-1170.0 




-1491.7 



Name 


Description 


Value 




specific membrane resistance 


30 krj-cm2 


R. 


specific axial resistance 


150 n-cm 


Cm 


specific membrane capacitance 


1.0 /xF/cm2 




Name 


Description 


Value 


r 


synaptic time constant 


0.1 ms 


le 


excitatory current amplitude 


50 pA 




Name 


Description 


Value 


^syn 


number of synapses 


1000 


^syn 


synaptic input rate 


5 spikes/s 


fill 


pairwiso injjut corrolatiou l)ot\voon colls 


{0.01.0.1.1.0} 




Name 


Description 


Value 


^ cond 


extracellular conductivity 


0.3 S/m 


F 






Name 


Description 


Value 


dtgim 


time resolution in simulation 


1/64 ms 


dtdata 


time resolution of data 


1.0 ms 


T 


simulation time 


either 10200 ms or 1200 ms 






(off-center) (^) 



(^) the somata of the L5 cells were adjusted 25 upwards compared to midpoint of L5 so that apical 
dendrites reach LI 

(^) the first 200 ms were discarded to avoid upstart effects 



Table S3. Parameters of the population model used for LFP simulations 
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